Coarsening Dynamics of Crystalline Thin Films 
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The formation of pyramid-like structures in thin-film growth on substrates with a quadratic sym- 
metry, e.g., {001} surfaces, is shown to exhibit anisotropic scaling as there exist two length scales 
with different time dependences. Analytical and numerical results indicate that for most realizations 
coarsening of mounds is described by an exponent n = l/(3v2). However, depending on material 
parameters, n may lie between (logarithmic coarsening) and 1/3. In contrast, growth on substrates 
with triangular symmetries ({111} surfaces) is dominated by a single length ~ t 1 ^ 3 . 

PACS numbers: 05.70.Ln, 64.60.My, 68.35.Fx, 68.55.-a 



43 

o 



CO 



I 

C 

O 

o 



On 



OO 
O 
OO 
On 



-a 
c 

o 
o 



X 



The formation of pyramids or mounds in growth pro- 
cesses such as molecular beam epitaxy has been observed 
in numerous experiments 0-^] and computer simulations 
[^|-^|. This three-dimensional growth mode is caused by 
step-edge barriers that hinder the diffusion of adatoms 
across step edges. This generates a surface diffusion cur- 
rent j that causes the initial high-symmetry orientation 
of the surface to become unstable |?],|2| . As a consequence 
the surface breaks up forming pyramid-like mounds, the 
side orientations of which correspond to stable zeros of 
the diffusion current ||J|^] . At this level this phenomenon 
is theoretically well understood. Conversely, the obser- 
vation that the evolving surface morphology coarsens as 
more and more material is deposited, is poorly under- 
stood. The similarity with ordering phenomena like 
spinodal decomposition has led to the assumption that 
correlation functions like the height-height correlation 
function 

C(r,t)=L- 2 J2(Hr + ^t)h( X ,t)) (1) 

X 

or its Fourier transform, the structure factor 



S(M) = (MM)M-M)) , 



(2) 



depend only on a single length scale, the pyramid size 
R(t), which in turn, because of the scale invariance of the 
correlation functions, must have a power-law dependence 
R(t) ~ t". Here, the surface height h(r,t) is measured 
in a comoving frame of reference such that the average 
height is zero, h(k,t) is the Fourier transform of h(r,t), 
and L 2 is the size of the substrate. In numerical simu- 
lations and several experiments the coarsening exponent 
was found to be close to n ~ 1/4. However, theories that 
attempt to calculate this exponent are so far not very 
convincing. 

Surprisingly, the initial assumption that there is only 
a single length scale and that the correlation function 
therefore follow simple scaling laws, was never really 
questioned. Only recently jl(^] it was found, that this 
kind of surface dynamics is in some respects differ- 
ent from ordering dynamics like spinodal decomposition: 
The dynamical exponent z that describes the relaxation 



of shape fluctuations of pyramids is not equal to the in- 
verse 1/n of the coarsening exponent indicating that the 
dynamical evolution in these two cases is governed by dif- 
ferent length scales. In this article it will be shown that 
the coarsening dynamics itself depend on at least two 
length scales with different time dependences. There- 
fore, the standard scaling assumption for the correlation 
functions is invalid. 

At not too high temperatures, so that desorption can 
be neglected, the evolution of the surface in molecular 
beam epitaxy is described by the continuity equation 



d t h = -AAh - V • j(m) + T) 



(3) 



where m = V/i is the slope of the surface profile and r\ 
represents shot noise due to fluctuations in the deposition 
flux. The form of the surface current j must be chosen 
such that it describes the above mentioned instability and 
leads to slope selection. The form j ~ m(l — m 2 ) has 
been used P,|Tl[, however, this expression has stable zeros 
for all slopes with |m| = 1 regardless of the direction of 
m. Such an azimuthal symmetry is unrealistic for crys- 
talline films. The simplest form that describes growth on 
substrates with a quadratic symmetry is a current with 
components 



jx = m x (l 



bm 2 ) , j y = m y (l - m 2 - bm 2 x ) (4) 



which leads to a buildup of pyramids with selected slopes 
(±1,±1)/(1 + b) for -1 < b < 1. We will concentrate on 
the case 6 = first. The relevance of the parameter b 
and triangular symmetries applicable to (111) substrates 
will be discussed further below. 

Coarsening in phase ordering dynamics, e.g., spinodal 
decomposition or Ostwald ripening, is a deterministic 
process, i.e., noise is irrelevant In the same spirit 

only the deterministic aspects of the dynamics described 
by Eq. (||) with n = will be studied here. Effects of the 
noise are discussed at the end of the article. 

The important differences between the problem stud- 
ied here and phase ordering dynamics described by the 
Cahn-Hilliard equation Jl3[ | becomes apparent, when the 
domain configurations are plotted as in Fig. [j]. A domain 
in this context is an area of constant slope corresponding 
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FIG. 1. a) Domain wall structure obtained from a numer- 
ical solution of Eqs. (^), (^); b, c) types of domain walls: 
b) pyramid edges, c) roof top. The greyscales correspond 
to average slopes in the following way: white = (+1, +1), 
blacks (-1, -1), light grey = (+1, -1), dark grey = (-1, +1). 



to one of the four values (±1,±1). The analogous case 
in phase ordering dynamics is described by a four-state 
clock model pi] . However, in that case domain walls 
do not have any particular orientation, whereas here do- 
main walls are intersections of planes of constant slopes 
and therefore form straight lines. Furthermore, there are 
two types of domain walls: Domain walls at which only 
one component of the slope changes are aligned along 
the x- and y-axes. These are the edges of the pyra- 
mids. Domain walls at which both components of the 
slope change run at 45° with respect to the principal 
axes. These latter domain walls form roof tops as illus- 
trated in Fig. 03. Domain walls in systems that phase 
order give rise to a power-law tail in the structure-factor 
called Porod's law [|||. On scales £ < x < R(t) the 
height function close to an edge of a pyramid running in 
the y-direction has the form h(x,y) = \x\ + y. Here, £ 
is the width of such a domain wall [£ = y/2 for (||)and 
6 = 0]. Hence, the slope m = (sign(x), 1) leads to a sin- 
gular contribution in the slope-slope correlation function 
@, Casing = —2\x\pnoo}, where /O(ioo) is the density 
of these domain walls jl7| . More generally, domain walls 
that make an angle ip with the y-axis yield a contribution 
S m = 4k~ 3 5(ip — Tp)p^{t) to the slope structure factor, 
where p^ (t) is again the density of the domain walls and 
tp the azimuthal angle such that k = (k cos ip, k simp). 
Consequently, the structure factor (j^) for crystalline films 
with a quadratic anisotropy has a highly anisotropic tail: 
for 1/R{t) < k < l/£ 



A 

s (M = fc5Eh°o>(*K 



ip — —v 



v=0 



-2p(iio) (t)S {v-\-\ v 



(5) 



where the additional factor 2 in front of the density 
P{no)W °f domain walls at 45° directions (roof tops) 
comes from the fact that both components of the slope 
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FIG. 2. Time dependence of the length scales l/fc(ioo)(t) 
and l/fyno) (t) as defined in Eqs. (^|), (Q). The average pyra- 
mid size R{t) measured as the first zero of the correlation 
function C(r,t) in the (110) directions is plotted as well. 

change at such a domain wall. Eq. (|J) is one of the 
central results of this article: It shows that the struc- 
ture factor depends on two different length scales 1 / /3(ioo) 
and l/jO(no), specifying the average separation of edges 
of pyramids and roof tops, respectively. Whereas for 
P(ioo) ^ P(no) the former is directly related to the aver- 
age pyramid size R(t), no obvious relation exists between 
P(iio) and R(t). In fact, there is no reason why l/p(ioo) 
and l/p(no) should have the same time- dependence or 
follow the same power law. Therefore, the structure fac- 
tor does not obey a scaling law S(k, t) = R 2 (t)s(kR(t)), 
not even in a modified anisotropic form fll8[ . Hence, all 
theoretical approaches that assume such a dependence 
on a single length scale, are unfounded. 

The result that the structure factor is nonzero only 
along the (100) and (110) directions as indicated by the 
delta functions in (|^) is a consequence of the assumption 
that the singular part of C(r, t) dominates the contri- 
butions to the tail of S(k,t). For directions other than 
these high symmetry directions no singularities exist and 
one should expect an exponential decay. Furthermore, 
although fluctuations of the domain wall directions are 
finite [^0| , such fluctuations will remove the delta sin- 
gularities and replace them by a sharply peaked func- 
tion with finite width. Nevertheless, the prediction that 
5(k, t) ~ fc -5 in directions normal to the domain walls, 
and that 5(k, t) decays faster for all other directions is 
confirmed in numerical simulations |h| and it should also 
be possible to observe this behaviour experimentally. 

Length scales that show the anisotropic scaling of the 
structure factor in a numerical solution of Eqs. (|^^) can 
be defined as inverses of moments of the structure factor 

fc <ioo) (*) = E fex k Mk x , 0), t) S((k x ,0), t) , (6) 

fc <no> (*) = E fc ^ k S((K k), t) /J2 k S((k, k), t) (7) 

in the two directions in question. Fig. || shows clearly 
that the scaling behaviour of S(k, t) depends on the di- 
rection of k: whereas l/fc(no)(i) is well described by a 
power law ~ t 1 ^ 3 , l/fc<ioo)(*) grows more slowly. In fact, 
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in a numerical simulation coarsening stops as soon as 
the average distance l/femo) between two roof tops is 
of the same order as the system size. At the same time 
the average pyramid size can be several orders of mag- 
nitude smaller. The coarsening of pyramids is enslaved 
to dynamics of roof tops: A configuration without any 
roof tops does not coarsen, i.e., it is a metastable state 
p9[ . Roof tops represent defects in such perfect pyramid 
lattices and coarsening proceeds by eliminating such de- 
fects. However, because of the more complicated scaling 
behaviour indicated in (|5|) the theory of Bray and Ruten- 
berg [^0| can no longer be used to calculate the growth 
law. In the limit where the average roof-top distance 
D r~j l/p(no) is much larger than 1//O(ioo) ~ R there ex- 
ists a simple geometric relation between the two length 
scales: Coarsening proceeds by elimination of finger-like 
domain configurations, two of which are indicated by ar- 
rows in Fig. |l|a. Such fingers are always terminated by 
roof tops of length ^/2R(t) at either end. If such a finger 
disappears within a time interval At, the roof top den- 
sity changes by A(l/D) = —2^/2R/L 2 , where L is some 
unit of length. In the same time the average density of 
pyramid edges changes by A(l/i?) = -2D/L 2 . It follows 
that £>A(l/£>) = y/2RA(l/R) or Z^i) 1 /^ _ R( t ). Us- 
ing the result D ~ t 1 / 3 (see below) one finds R ~ t n with 
n = l/(3y/2) ~ 0.2357, hardly distinguishable from 1/4 
in a numerical simulation. 

During the temporal evolution of the surface morphol- 
ogy shown in Fig. |l| many more pyramids than roof tops 
are formed. The reason lies in the value 6 = that 
was used in that simulation. The parameter b describes 
how the rr-component of the current depends on the y- 
component of the slope. Such transverse currents are 
typical for effects like edge diffusion, the importance of 
which for growth morphologies has been emphasized by 
Bartelt and Evans pi]] and was recently confirmed in an 
experiment p^] . It is possible to assign a "surface ten- 
sion" a to domain walls by integrating over the domain 
wall profile, 



dr i 



ch 



drj_ 



dr± 



(8) 



where the coordinate r± runs perpendicular to the do- 
main wall and the domains are assumed to extend to 
infinity on both sides of the wall. Using the steady-state 
solutions for the domain wall profiles that correspond 
to Eqs. (§), (|), it is found that a pc = 4/[3y / 2(l + b)} 
for pyramid edges and er r t = 4/[3\/l + b] for roof tops. 
Therefore, for b — the formation of roof tops is sup- 
pressed. To understand the coarsening mechanism it is 
instructive to study perturbations to a perfect lattice 
of pyramids (see Fig. ||a-e). The movement of domain 
walls necessarily involves the nucleation of roof tops, see 
Fig. pp, c. Such a perturbation does not cost anything, 
if v^o"rt = Cpc or, if b = —3/4. For this value of the 
parameter b edge-like domain walls can transform with- 
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FIG. 3. Domain wall motion for films with a quadratic 
anisotropy: a-e) schematics of domain wall motions through 
the introduction of roof tops, f) coarsening dynamics for 
the special degenerate case b = —3/4. Data for triangular 
anisotropies (triangles) are included as well. The straight lines 
correspond to power laws ~ t 1 ' 3 . 



out cost into roof-top-like domain walls and coarsening 
is no longer enslaved to the dynamics of the roof tops. 
This is convincingly confirmed by numerical solutions of 
Eq. (0) with a current (|4|) and b = —3/4. As seen in 
Fig. pp coarsening is much faster as all length scales in- 
crease as i 1 / 3 , i.e., the coarsening of the pyramid size R(t) 
depends on the parameter b that in turn is determined by 
microscopic processes like edge diffusion. 

The time scale for the elimination of roof tops is set by 
the time it takes to get from configuration ||c to configura- 
tion |e. This is a simple diffusion process as the changes 
indicated in ||d do not change the lengths of the domain 
walls. This behavour and the time dependence ~ t 1 / 3 
suggest that an approach similar to the Lifshitz-Slyozov- 
Wagner theory ^3| can be used here. Assuming that the 
magnitude of the slope is slightly less than the asymptotic 



value, i.e. 



(r,t)| 



m 



(r, t) with e X)3/ <C mo 



it can be shown that there exists a diffusion current from 
the small domains to the larger domains that leads to 
a coarsening law ~ i 1 / 3 . However, the same theory also 
predicts that the average deviation (e(t)} should decrease 
as t -1 / 3 . This latter behaviour is not confirmed in nu- 
merical solutions of Eq. (p|) for triangular symmetries or 
quadratic symmetries with b = —3/4 [ ^9[ . This problem 
requires further investigation. 

In the nondegenerate case with b > —3/4 the configu- 
rations ||b, c are highly suppressed as they correspond to 
activated processes. Hence, coarsening proceeds by elim- 
inating roof tops that were nucleated in the initial stages 
of the instability. Consequently, the roof-top distance 
becomes much larger than the pyramid size R. 

Although coarsening exponents close to 1/3 were found 
in a numerical simulation the situation that the sur- 
face tensions obey the relation v2cr r t = a pc exactly is 
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unlikely to be realized in nature. The fact that many ex- 
periments measured a coarsening exponent close to 1/4 
instead indicates that typically the pyramid coarsening 
is enslaved to the dynamics of the roof tops in which case 
the theoretical value of the exponent is l/(3\/2) as ex- 
plained above. However, this is no longer correct for thin 
films with a triangular lattice anisotropy, i.e., for (111) 
substrates. In that case there is only one type of domain 
wall; thus one expects fast coarsening in any case. The 
appropriate form of the surface current in that case was 
given in Ref. Jl4|]. Results of a numerical solutions of 
Eq. (|^) are included in Fig. |^f (triangles). They indeed 
show that the pyramid size R(t) ~ i 1 / 3 . The author is 
aware of a single experiment that studied coarsening on 
a (111) substrate [Q. The results of that experiment 
are in perfect agreement with the theory presented here. 
Nevertheless, it must be emphasized that even in those 
cases where coarsening is ~ t 1 / 3 the correlation functions 
( 1 ) and (^) do not obey a simple scaling law. The result 
(5) that the structure factor scales differently in direc- 
tions perpendicular to the domain walls than in all other 
directions is still valid. 

Most of the results derived in this article do not de- 
pend on the specific form of the surface current, e.g., 
the anisotropic tail of the structure factor (||) and the 
growth exponents follow from the constraints imposed 
by the crystalline anisotropics of the growing film, which 
in turn severely restrict the possible domain configura- 
tions. Only the results for <7po and <r r t are specific to the 
form of the surface current (0) . Because of the existence 
of several length scales and the fact that the correlation 
functions do not obey a simple scaling law, the method by 
which the characteristic length scales are measured be- 
comes important: Whereas the average roof-top distance 
could be separated from the average pyramid size using 
moments of the structure factor in different directions 
(||), (0), this cannot be accomplished in real space, e.g., 
the height-height correlation function (Q) in the (100) di- 
rections changes its functional form with time and does 
not permit the determination of any length scale. 

As mentioned in the introduction, most experiments 
measured a coarsening exponent for the pyramid size on 
substrates with quadratic symmetries of n < 1/4 indicat- 
ing that real systems are nondegenerate, i.e., they corre- 
spond to b > —3/4 in Eq. (Q). In these cases the roof-top 
distance D becomes much larger than the pyramid size 
R and R ~ t 1 /^) wne reas D ~ t 1/3 . However, if the 
late stages of the coarsening process are studied on scales 
smaller than D, the pyramid size coarsens only logarith- 
mically as is typical for activated processes and as has in 
fact been predicted for these kind of systems ]25| ]. The 
measurement of the coarsening behaviour in these cases 
is severely complicated by a wide cross-over regime |i~9| ] 
where almost any value between n = and n < 1/4 
can be observed. The interpretation is further compli- 
cated when effects of the noise r\ ^ is included: For 
some time such stochastic fluctuations provide a mech- 



anism to overcome the activation barriers that suppress 
the formation of roof tops. However, the activation bar- 
riers are proportional to the pyramid size. Hence, in the 
asymptotic regime, where the pyramid size is sufficiently 
large, stochastic effects become unimportant. This is in 
agreement with arguments presented by Shore et al. p5j . 
For deposition on (111) substrates noise is clearly irrele- 
vant. 

The author thanks M. Plischkc for many helpful con- 
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of Canada. 
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